Nucleon mass from a covariant three-quark Faddeev equation 
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We report the first study of the nucleon where the full Poincare-covariant structure of the three-quark 
amplitude is implemented in the Faddeev equation. We employ an interaction kernel which is consistent 
with contemporary studies of meson properties and aspects of chiral symmetry and its dynamical breaking, 
thus yielding a comprehensive approach to hadron physics. The resulting current-mass evolution of the 
nucleon mass compares well with lattice data and deviates only by ~ 5% from the quark-diquark result 
obtained in previous studies. 
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Experiments, and hereby especially electroweak probes 
at all energy scales, have provided detailed information 
about the structure of the nucleon. Nevertheless, under- 
standing the nucleon's structure in terms of quarks and glu- 
ons, the elementary degrees of freedom of Quantum Chro- 
modynamics (QCD), has remained a challenge in theoreti- 
cal hadron physics. 

Starting with the original work of Faddeev [Q]] a formal- 
ism has been developed to treat a relativistic three -body 
problem In its covariant form, the Faddeev equation 

is the three -body analogue of the two-body Bethe-Salpeter 
equation (BSE) |5fl. In the case of the nucleon, its solution 
is a covariant three-quark amplitude whose relativistic spin 
structure has been explored in Q] and in the light-front 
formalism in is- Hlh . A complete classification according 
to the Lorentz group and the permutation group §> 3 was de- 
rived in lfl2ll in terms of covariant three- spinors. 

The formalism of QCD's Dyson-Schwinger equations 
(for recent reviews, see e.g. Ill3lll4ll ) provides a way to em- 
bed the covariant three-quark Faddeev equation in a consis- 
tent quantum-field theoretical setup. The dynamical ingre- 
dients in the equation (the dressed quark propagator and 
the three-quark kernel) can then be treated in perfect cor- 
respondence with studies of quark and meson properties as 
well as related aspects of QCD. 

The biggest obstacle on the way to a direct numerical 
solution of the three -body bound-state equation is its com- 
plexity. Simplifications employed in the past implemented 
perturbative quark propagators lfl5l flr^l . together with a 
three-body spectator approximation ltl7ll . or in a Salpeter- 
equation setup with instantaneous forces 0]. The corre- 
sponding equation of a scalar three-particle system with 
scalar two-body exchange was recently investigated and 
compared to the light-front approach Il8ll . Another kind of 
simplification can be achieved by considering diquark cor- 
relations (see e.g. lfl9h for an overview). While maintaining 
full Poincare covariance, the quark-diquark model traces 
the nucleon's binding to colored scalar- and axialvector 
diquarks, thereby simplifying the Faddeev equation to a 
quark-diquark BSE. This strategy has been applied to study 
nucleon and A properties Il2fj| - l23ll . 



Here we report the first fully Poincare-covariant com- 
putation of the nucleon's Faddeev amplitude beyond the 
quark-diquark approximation. The numerical solution of 
the Faddeev equation is performed after truncating the in- 
teraction kernel to a ladder dressed-gluon exchange be- 
tween any two quarks, thereby enabling a direct compar- 
ison with corresponding meson studies as well as earlier 
investigations of baryons in the quark-diquark model. 

In QCD baryons appear as poles in the three-quark scat- 
tering matrix. This allows one to derive a relativistic three- 
body bound-state equation: 



\£ = K (3) * 



Here, is the bound-state amplitude denned on the 
baryon's mass shell. The three -body kernel comprises 
a three-quark irreducible contribution and the sum of per- 
muted two-quark kernels whose quark- an tiquark analogues 
appear in a meson BSE, and the superscript a denotes the 
respective accompanying spectator quark. 

The observation of a strong attraction in the SU (3)c an- 
titriplet qq channel has been the guiding idea for the quark- 
diquark model, namely that quark-quark correlations pro- 
vide important binding structure in baryons. This moti- 
vates the omission of the three-body irreducible contribu- 
tion from the full three-quark kernel. The resulting covari- 
ant Faddeev equation includes a sum of permuted qq ker- 
nels (cf. Fig.[T|: 
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where denotes the renormalization-group invariant 

products of a qq kernel and two dressed quark propagators: 

act 

{a, b, c} is an even permutation of {1, 2, 3} and linked to 
the respective Dirac index pairs. 




The spin-momentum part of the full Poincare-covariant 
nucleon amplitude ^^(j?, q, P) is a spin- 1 ^ four-point 
function with positive parity and positive energy: it car- 
ries three spinor indices {a, j3, 7} for the involved valence 
quarks and one index 5 for the spin-1 /2 nucleon. The am- 
plitude depends on the total momentum P and two relative 
Jacobi momenta p and q, where P 2 = —M 2 is fixed. It 
can be decomposed into 64 Dirac structures: 



64 



^ aPl s{p, q,P) = J2 h r^ jS {P: 9, p )> 



(4) 



k=l 



where the amplitude dressing functions f k depend on the 
five Lorentz-invariant combinations 



P 



z 



Pr-qr, Zi=p-P, z 2 = q-P. (5) 



Here, a hat denotes a normalized 4- vector and pt^ = Tf VanV 



a transverse projection with Tp 1 " = 5 fJll/ — P^P V . 

A general spinor four-point function which depends on 
3 independent momenta involves 128 independent compo- 
nents of positive parity. An orthogonal basis {r k } for the 
64-dimensional subspace of a positive-parity and positive- 
energy nucleon is given by the set 



Dr 



1 (8)1 
7 5 ® T 



(T i ®r j )(A r l5 C®A+), (6) 



7 4 7 2 is the charge-conjugation matrix, r 



where C 

± refers to the positive- and negative-energy projectors 
A ± (P) = (1 ± P)/% an d the tensor product is under- 



0-yS 



The relative-momentum 



stood as (A(&B) a 
dependence of the basis elements is carried by the T it 
i = 1,2, 3, 4, defined by 



Ti(p,q,P) = {l, \ \p T Atl Pt, 4t} ■ 



(7) 



The momenta {pr, q f , P} were conveniently chosen to be 
orthonormal with respect to the Euclidean metric via 

p£:=r£V, qr.= T^T» x q x =T^q». (8) 

A partial-wave decomposition leads to linear combina- 
tions of the {Slj, Py} as eigenstates of quark-spin and or- 
bital angular momentum operators S 2 and L 2 in the nu- 
cleon rest frame. The 64 basis covariants (32 each for total 
quark spin s = 1/2 and s = 3/2, respectively) can be ar- 
ranged into sets of 8 s-waves (I = 0), 36 p-waves (I = 1), 



and 20 (f-waves (I = 2). For instance, the dominant contri- 
butions to the Faddeev amplitude are given by the s-waves 



7sC ® A+ = ]T 



(9) 

^ T C ® 7^7 5 A + = £ (r S^ 2 + + P\ A ) , 

with 7^ = Tp V ^ y '. In the quark-diquark model, these cor- 
respond to scalar-scalar and axialvector-axialvector combi- 
nations of diquark and quark-diquark amplitudes for either 
of the three diagrams appearing in the Faddeev equation. 

The basis elements can be expressed in terms of 
quark three-spinors frequently used in the literature, e.g. 
Ref. Ull. In this context the elements = A + 7 5 C(g)A + 
and := A+7^C ® 7t7sA + read 



-S+ W 



(Wu l + u l W)W 



(10) 



where the If 7 (P) are eigenspinors of A + and therefore sat- 
isfy the free Dirac equation for a spin- 1 ^ particle. 

The Pauli principle requires the Faddeev amplitude to be 
antisymmetric under exchange of any two quarks. The Fad- 
deev kernel Km is invariant under the permutation group 
§ 3 . The eigenstates of the Faddeev kernel can hence be 
arranged into irreducible S 3 multiplets 



A, 



M A 
Ms 



(ID 



of which the first two (totally symmetric or antisymmetric) 
solutions are unphysical while the mixed-symmetry dou- 
blet constitutes the Dirac part of the nucleon amplitude. 
Taking into account the flavor and color structure, the full 
Dirac-flavor-color amplitude reads 



^{p,q,P) = [^ Mj Jm a +^> 



M s 1 M. 



1 £AB£ 



(12) 



where T Ma , T Ms denote the isospin-1/2 flavor tensors 
for proton and neutron and Eabc the antisymmetric color- 
singlet wave function. A flavor-dependent kernel in the 
Faddeev equation will mix ^ m a anc ^ ^ m s whose dom- 
inant contributions are given by and Af v respectively. 
Similarly to the analogous case of a diquark amplitude, the 
symmetry does however not reduce the number of Dirac 
covariants since the dressing functions f k transform under 
the permutation group as well. 
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To proceed with the numerical solution of the Faddeev 
equation, we need to specify the quark-quark kernel fC and 
the dressed quark propagator S(p) which appear in Eq. ([3]). 
This is achieved via the axial- vector Ward-Takahashi iden- 
tity which encodes the properties of chiral symmetry in 
connection with QCD. Its satisfaction by the interaction 
kernels in related equations guarantees the correct imple- 
mentation of chiral symmetry and its dynamical breaking, 
leading e.g. to a generalized Gell-Mann-Oakes-Renner re- 
lation valid for all pseudoscalar mesons and all current- 
quark masses 1124 12511 . In particular the pion, being the 
Goldstone boson related to dynamical chiral symmetry 
breaking, becomes massless in the chiral limit, independent 
of the details of the interaction. Specifically, we describe 
the qq kernel by a ladder dressed-gluon exchange: 



fcaa'0j3'{k) — Z. 



2 47ra(fe 2 ) 



A: 2 



TTtI'T^ (13) 



which must also appear in the corresponding quark Dyson- 
Schwinger equation whose solution defines the renormal- 
ized dressed quark propagator: 



Z 2 (ip + 



(14) 

The bare quark mass m enters as an input, and the gluon 
momentum is k = p — q. The inherent color structure of 
the kernel leads to prefactors 2/3 and 4/3 for the integrals 
in Eqs. d23> an d ( fl4l . respectively. 

Eqs. ( TT31 - IT4T ) define the rainbow-ladder (RL) truncation 
which has been extensively used in Dyson-Schwinger 
equation studies of mesons and baryons in the quark- 
diquark model, e. g. ll26l l27ll and references therein. The 
non-perturbative dressing of the gluon propagator and the 
quark-gluon vertex are absorbed into an effective coupling 
a{k 2 ) for which we adopt the ansatz [28, 29h 



a(k 2 



A 2 



+ auv(k 2 



(15) 



The second term reproduces the logarithmic decrease of 
QCD's perturbative running coupling and vanishes at k 2 = 
0. The first term is parametrized by an infrared scale A and 
a dimensionless parameter r]. It yields the non-perturbative 
enhancement at small and intermediate gluon momenta 
necessary to generate dynamical chiral symmetry break- 
ing and hence a constituent-quark mass scale. ({A, 77} 
and the infrared parameters used in [29] are related by 
C = (A/A t ) 3 and u = r^A/At, with A t = 1 GeV.) 

Beyond the present truncation, corrections arise from 
pseudoscalar meson-cloud contributions which provide a 
substantial attractive contribution to the 'quark core' of 
dynamically generated hadron observables in the chiral 
regime and vanish with increasing current-quark mass, but 
also from non-resonant contributions due to the infrared 
structure of the quark-gluon vertex. To anticipate correc- 
tions we exploit the freedom in adjusting the input scale A. 





Q-DQ[23J 


Faddeev (Mj\) 


Faddeev (A^) 


Setup A 


0.94 


0.99 


0.97 


Setup B 


1.26(2) 


1.33(2) 


1.31(2) 



TABLE I. Nucleon masses obtained from the Faddeev equation 
in setups A and B and compared to the quark-diquark result. The 
77 dependence is indicated for setup B in parentheses. 



We adopt two different choices established in the literature 
in the context of tt and p properties II29I1 : 

Setup A is determined by a fixed scale A = 0.72 GeV, 
chosen in [28] to reproduce the experimental pion decay 
constant and the phenomenological quark condensate. Cor- 
responding results are therefore aimed in principle at a 
comparison to experimental data for meson and baryon 



properties (see 11231 12611 and references therein). Setup B 
defines a current-mass dependent scale which is deliber- 
ately inflated close to the chiral limit, where A ~ 1 GeV 
ll29ll . It is meant to describe a hadronic quark core which 
must subsequently be dressed by pion-cloud effects and 
other corrections. As a result, tt, p, N and A observables 
are consistently overestimated, but (with the exception of 
the A-baryon) compatible with quark-core estimates from 
quark models and chiral perturbation theory (for a detailed 
discussion, see i22l. I23L, l29lo . Irrespective of the choice of 
A, hadronic ground-state properties have turned out to be 
insensitive to the value of 7/ in a certain range [26l 281. 
Consequently, with Eqs. ([T3l4X5T l and A, the input of the 
Faddeev equation is completely specified with all parame- 
ters already fixed to meson properties. 

Since the dressed-gluon exchange kernel is flavor- 
independent and we consider only equal quark masses, 
the equations for the Dirac amplitudes ^ m a anc ^ ^ m s m 
Eq. ([Lit decouple because of the orthogonality of the two 
flavor tensors T m a anc ^ ^ Ms- H ence one obtains two de- 
generate solutions of the Faddeev equation, where by virtue 
of the iterative solution method the symmetry of the start 
function determines the symmetry of the resulting ampli- 
tude. The massive computational demand in solving the 
equation primarily comes from the five Lorentz-invariant 
momentum combinations of Eq. (0 upon which the ampli- 
tudes depend. In analogy to the separability assumption of 
the nucleon amplitude in the quark-diquark model we omit 
the dependence on the angular variable z a = p~r ■ qr but 
solve for all 64 dressing functions fk(p 2 , q 2 , 0, Zi, z 2 ). 

The resulting nucleon masses at the physical pion mass 
in both setups A and B are presented in Table U The dif- 
ference of ~ 2% between the Ma and Ms solutions is 
presumably an artifact associated with the omission of the 
angle z . For either solution typically only a small number 
of covariants are relevant which are predominantly s-wave 
with a small p-wave admixture. The angular dependence 
in the variable z 2 is small compared to Z\ in analogy to the 
quark-diquark model, where the dependence on the angle 
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FIG. 2. (color online) Evolution with of m p and Mjv com- 
pared to lattice data; see (23ll for references. The quark-diquark 
model result for M N is plotted for comparison. Dashed and 
dashed-dotted lines correspond to setup A; the solid line for m p 
and the bands for Mjy (mixed-antisymmetric solution) are the re- 
sults of setup B, where the variation with -q is explicitly taken into 
account. Dots denote the experimental values. 



between the relative and total momentum of the two quarks 
in a diquark amplitude is weak. 

The evolution of M N and the p-meson mass from the 
BSE vs. m% is plotted in Fig. (0 and compared to lattice 
results. The findings are qualitatively similar to those for 
m p : setup A, where the coupling strength is adjusted to 
the experimental value of f w , agrees with the lattice data, 
which is reasonable in light of a recent study of corrections 
beyond RL truncation for the p- meson [30]. Setup B pro- 
vides a description of a quark core which overestimates the 
experimental values while it approaches the lattice results 
at larger quark masses. 

A comparison to the consistently obtained quark-diquark 
model result exhibits a discrepancy of only ~ 5%. This 
surprising and reassuring result indicates that a descrip- 
tion of the nucleon as a superposition of scalar and axial- 
vector diquark correlations that interact with the remain- 
ing quark provides a close approximation to the consistent 
three-quark nucleon amplitude. 

We have provided the first fully Poincare-covariant 
three-quark solution of the nucleon's Faddeev equation. 
The present study contains the first numerical results for 
the nucleon mass in this approach. Due to the consider- 
able computational efforts involved, more results and an 
in-depth investigation with regard to the complete set of 
invariant variables will follow in subsequent publications. 
Future extensions of the present work will include an anal- 
ogous investigation of the A-baryon, more sophisticated 
interaction kernels, e.g. in view of pionic corrections, and 
ultimately a comprehensive study of baryon resonances. 
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